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Abstract A hierarchical interval subdivision is shown 
to lead to a p-adic encoding of image data. This allows 
in the case of the relative pose problem in computer 
vision and photogrammetry to derive equations having 
2-adic numbers as coefhents, and to use the Hensel lift- 
ing method to their solution. This method is applied to 
the linear and non- linear equations coming from eight, 
seven or five point correspondences. An inherent prop- 
erty of the method is its robustness. 

Keywords Relative pose • p-adic numbers • Essential 
matrix ■ n-point method 



1 Introduction 

The issue of estimating camera motion from two views 
is classical by now, and methods from projective and al- 
gebraic geometry towards its solution were employed at 
an early stage (e.g. p~5|). The beginning of this present 
century witnesses the application of sophisticated meth- 
ods from computational commutative algebra in order 
to rephrase the equations into a form from which so- 
lutions can be obtained with relative ease. The rela- 
tionship between the views is established by finding 
correspondences between point pairs taken from both 
images. The fundamental matrix faithfully encodes the 
geometric relationship between the two images. For nor- 
malised cameras, the fundamental matrix coincides with 
the essential matrix. In general, the two matrices are 
related through the camera calibration. Hence, if the 
calibration is known, it is sufficient to estimate the es- 
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sential matrix in order to solve the relative pose prob- 
lem. From a conceptional as well as a computational 
point of view, it makes sense to use only few correspon- 
dences of image points in order to estimate the essential 
matrix E. And different samples of n correspondences 
lead to a set of candidate essential matrices from which 
an optimal choice can be made. 

Since E — (e^) is a projective 3x3 matrix, i.e. 
only determined up to a scalar factor, the number of 
parameters to be found is at most 8. As each pair of 
corresponding image points leads to a linear equation 
in the e^, it suffices to take n — 8, as has been showed in 
pTj . However, this linear method does not take into ac- 
count further constraints on E. So, a 7-point algorithm 
was developped [6l[9] (cf. also \\7 j for an overview). 
The 5-point algorithm of |I3j uses the minimal number 
of point correspondences necessary for determining E. 
The non-linear constraints lead to homogeneous equa- 
tions of degree three in four unknowns which are trans- 
formed through an elimination process into a univariate 
equation of degree 10. Hence the number of complex so- 
lutions is not more than 10. The elimination process can 
be simplified [TU] , or replaced by Grobner basis methods 
|14j . There is a lot of existing work towards optimising 
the performance of this method in the computational 
algebraic geometry community. Stewenius et al. com- 
pare the performance of elimination and Grobner basis 
methods [16J. 

Correspondences often suffer a geometric perturba- 
tion, i.e. the correct point u' in image I 1 corresponding 
to point u in image I is mostly found up to a small er- 
ror vector e such that u! + e instead of u' is matched to 
u. In this article, a hierarchical method based on inter- 
val subdivision is developped to handle this problem. A 
natural way of encoding hierarchies is provided by the 
p-adic numbers, where p is a fixed prime number. In 



2 



our case, we can use p = 2 and represent image points 
by pairs of binary expansions 

oo 

a = >J fli 2 * 

with coefficients m equal to or f. These expansions 
can be infinite, theoretically. Practically, the finiteness 
of resolution means approximation through truncation. 
The framework for this method is p-adic geometry which 
has been applied in video segmentation and data anal- 
ysis pp2p] . 

In the context of relative pose problems, we will use 
the so-called lifting method for solving the equations. 
This is provided by Hensel's lemma which says that 
under certain conditions a solution of a given equation 
modulo p can be expanded to a p-adic solution. Its proof 
is constructive. In fact, this is the p-adic analogue of a 
Newton iteration. Applied to the equations of the n- 
point problems for n = 5, 7, 8, we obtain the result that 
for many choices of point correspondences, a 2-adic so- 
lution can be constructed. The 2-adic essential matrix 
then allows to hierarchically approximate the rigid mo- 
tion by truncation. A side effect of the lifting method 
is its high robustness to geometric perturbations. The 
encoding method ensures further that the number of 
Newton iterations required is proportional to the order 
of resolution desired. 

The article is structured as follows. The next two 
sections are a brief introduction to p-adic numbers, and 
to Hensel's lemma in a multivariate formulation. This 
is followed by a section on 2-adic encoding of image pix- 
els for a p-adic camera model. The last section applies 
Hensel's lemma to the 2-adically defined equations for 
the eight-, seven-, and five-point problems. 

2 p-adic numbers 



sums 

aN = ^ aipl 

converge to a: 

\a - a N \ p = J J +1 — >0 for N — > 00. 

The domain of all p-adic numbers is denoted by Q p and 
densly contains the rational numbers Q with respect to 
this metric. Those p-adic numbers a with v a > are the 
p-adic integers, denoted as Z p . This domain contains 
densely the usual integers Z. An equivalent description 
of p-adic integers is given by 

Z p = {a e Q p I \a\ p < l}. 

Approximation of p-adic integers by their partial 
sums, simply termed "truncation" , has an algebraic for- 
mulation. In its simplest form, a p-adic integer a can be 
given by its coefficient a G {0, . . . ,p — I}, and another 
p-adic integer b having the same coefficient bo — ao ap- 
proximates a up to that order of magnitude. This is 
the case if and only if a — b is divisible by p. Hence, we 
arrive at 

Zp/pZp ^ Z/pZ = Fp, 

where the latter is the finite field with p elements. Sim- 
ilarly, the algebraic formulation of approximation up to 
terms of higher order is given by 

Zp/p^Zp = Z/p w Z =: Z p n. 

In other words, congruences modulo p N yield finite ap- 
proximations of p-adic numbers by their partial sums up 
to terms of order N. It is this property which makes p- 
adic numbers very suitable for hierarchically organised 
data. In later sections, we will see how this algebraic 
formulation can be used in stereo vision. A standard 
reference for p-adic numbers is [5] . 



The p-adic numbers were first described by Hensel in 
[5] . They are expansions of the form 

a = arf (I) 

into powers of a fixed prime number p and coefficients 
cij G {0, . . . , p — f }. If there are only finitely many terms 
in fl)). this defines a rational number. Any natural num- 
ber has a finite p-adic expansion (TTJ) with i > 0. The 
important observation is that expansions with infinitely 
many positive powers of p have a meaning. Namely, by 
defining 

where v a G Z is the smallest exponent occurring in 
expansion ([1]) , one obtains a metric for which the partial 



3 Hensel's lemma 

An important method in p-adic analysis is the so-called 
"lifting" of zeros of polynomials from F p to Z p . This 
uses the fact that 

Zp/pZp ^ Z/pZ = Fp 

by reducing a given equation f(X) = over Z p mod- 
ulo p to an equation / mod p over F p . Hensel's lemma 
then gives a criterion when a solution x\ of / mod p 
leads to a solution of /. This solution £, if it exists, is 
constructed by an iterated process. Namely, first a so- 
lution X2 of / mod p 2 is constructed from x\ , and from 
this a solution X3 of / mod p 3 etc. Each step of this 
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iteration yields an approximation to the true solution 
in Z p in the sense that 

f(xi) = mod p\ 

which in terms of the p-adic norm translates to 

1/(^)1 < - r 

In other words, the sequence f(x{) converges p-adically 
to the value /(£) = 0. The construction process guaran- 
tees further that Xi also converges to the p-adic solution 

t 

We now state a multivariate Hensel's lemma, but 
not in the most general form. For this, r L\X\ 1 . . . ,X n ] 
denotes the polynomial ring in n variables with integer 
coefficients. 

Theorem 1 (Multivariate Hensel's lemma) Let 

/ = (/ 1) ..,/ m )eZ[x 1) ...,x„r 

with m < n, and let k > 2. Suppose that the vector 
x = [x\, . . . , x n ) S 27 l is a solution of the congruences 

f 1 (X 1 ,...,X n )=0 mod/" 1 

f m (X 1 ,...,X n )=0 mod/" 1 
and that the matrix 




is modulo p ofrankm. Then there existt = [t\, . . . ,t n ) 6 
{0, . . . ,p— I}™ such that 

f{x + p k -H) = mod/. 

In particular, the equations f{X) — have a solution 
in Z p . 

Proof Consider the linear part in the Taylor expansion 
of / in x: 

f(X) = f(x) + Df(x) ■ (X — x) + terms of higher order, 

where X = (Xi, . . . ,X n ) denotes the vector of vari- 
ables. By assumption, it holds true that 

f(x) = -a 

for some vector a £ Z m . Due to the rank condition, the 
system of congruences 

a + D f (x) - t = modp (2) 
has a solution vector t 6 {0, . . . ,p— 1}™. Hence, 
f(x+p k -H) = f(x)+D f (x)-p k -H 

= p k -\a + D f {x) -t) =0 mod/. 
Iteration proves the last assertion. 



Usually, Hensel's lemma is stated in the case of a sin- 
gle univariate polynomial / having modulo p 1 a zero 
x. The rank condition translates to f'(x) ^ mod p, 
i.e. a; is a simple zero of / modulo p. The proof is a p- 
adic analogue of a Newton iteration. From ([2]) it follows 
for m = n that the rank condition implies uniqueness 
of the lift. The univariate case can be found e.g. in [5]. 

4 p-adic projective cameras 

4.1 p-adic encoding of image pixels 

Assume that a rectangular 2D image / is given by M x 
N pixels and that m, n are minimal such that M < 2 m , 
N < 2™. We will use a binary encoding of a given image 
point with pixel coordinates 

(x,y) e{0,...,M-l}x {0,...,N-1}. 

This means that x and y will be represented by a pair 
of binary numbers obtained in a hierarchical manner. 
Namley, consider one coordinate at a time, say x. It 
can be arrived at by a sequence of iterated subdivi- 
sions of the interval < • ■ • < 2™ — 1 into intervals of 
equal length. After m iterations each interval obtained 
contains precisely one pixel x-coordinate. This means 
that a given value of x lies in a uniquely determined 
nested sequence of intervals produced by this subdivi- 
sion process. The intervals form a rooted binary tree 
whose leaves correspond uniquely to the x-coordinates 
of pixels. By assigning for a given interval its left half 
the value 0, and 1 for its right half, we obtain the binary 
representation 

m — 1 

r(x) = a v r > c^e{0,l} 

i/=0 

by traversing the path from root down to x and picking 
up the zeros and ones along the way. 

We will interpret the natural number r(x) as a 2- 
adic integer: r{x) € 1^2- Each partial sum 

e 

i/=0 

is a 2-adic approximation of r{x). Its error is bounded 
by the 2-adic distance 

\n - r(x)\ 2 . 

In fact, r(x) itself is a 2-adic approximation to some 
imaginary "pixel" £ of infinite precision. The number 
m is given by the degree of resolution with which £ is 
viewed on the given image. And an infinitely precise 
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pixel would have a 2-adic expansion with possibly in- 
finitely many coefficients equal to 1. 

We proceed similarly for the y-coordinate, and ob- 
tain an encoding 

c 2 : I ^Z 2 xZ 2 , (x,y) ^ (r(x),s(y)), 

where s is the binary encoding of the y-coordinate. 

Any p-adic encoding of a number x bears the prob- 
lem that its p-adic approximations re are determined by 
arithmetic properties of x, ignoring their (euclidean) ge- 
ometric properties. This means that although xg could 
be p-adically very close to x, it could be at large dis- 
tance for the euclidean norm. Hence, arbitrary p-adic 
encodings of pixel data would be quite unsuitable for 
many applications in computer vision. However, the bi- 
nary encoding through interval splitting does not suffer 
this disadvantage. Namely, the map c 2 respects the im- 
age geometry in the following sense: 

Proposition 1 The map 

oo oo 

yields an inclusion of Z 2 into the real interval [0,2). It 
has the property 

\r-s\ 2 <2- i ^\i 2 (r)-i 2 (s)\ oc <2- e 

for all eem. 

Proof Any sum of negative powers of 2 can be ma- 
jorised by a (possibly shifted) geometric sum. Hence, 
such a sum converges for the real metric. This implies 
that the map l 2 is well-defined for arbitrary 2-adic num- 
bers. 

The inclusion property is clear, at least for finite 2- 
adic expansions. However, the map t 2 is injective on all 
of 1i2- Otherwise, assume that 

L 2 (a) = L 2 (h) 

for some a ^ b. Then a = ^2 an d b = b%2 1 must 
differ in some coefficient: 

Let n be minimal such that this occurs. Then it follows 
that 

|t 2 (a)-i 2 (6)| 0O > pr, 

because the inequality holds true for all partial sums 
with at least N terms. This is a contradiction. Hence, 
l 2 is injective. 

ssume |r — s| 2 < 2 e . This means that the 2-adic ex- 
pansions of r and s have the first £+1 terms in common. 
The last assertion follows from this. 

As a consequence, we obtain an embedding of the 
image / into [0, 2) 2 via composition of c 2 and i 2 x i 2 . 



4.2 p-adic camera model 

Recall that a projective camera is a map between pro- 
jective spaces 

k : P 3 -> P 2 

given by a 3 x 4-matrix of rank 3. Usually, cameras are 
modelled as being defined over the real numbers, i.e. 
they are given by real 3 x 4-matrices. However, when 
dealing with data coming from 3D-objects from the real 
world, cameras are usually approximated by rational 
matrices. In this way, we arrive at rational cameras: 

k: P 3 (Q) ^P 2 (Q) 

as maps between the rational points of projective spaces. 
For such cameras, the real completion is only one choice 
of many. Hence, p-adic cameras 

k: P 3 (Q p ) -» P 2 (Q P ) 

could also be considered in order to use methods from 
p-adic geometry for camera computations. 

In fact, since finite resolution prevents stereoscopic 
vision at arbitrarily large distance, it can become possi- 
ble to shift coordinate systems in such a way that cam- 
eras are described by matrices whose entries are natural 
numbers, at least approximately. In the following sub- 
section, we will show how this said approximation can 
be done for the 2-adic norm. Precisely the hierarchical 
method used in the previous subsection in 2D extends 
in a natural way to 3D in order to arrive at 2-adic cam- 
era model 

k: P 3 (Z 2 ) ^P 2 (Z 2 ). 

4.3 How to interpret a 2-adic essential matrix 

Assume that through correspondences between some 
points u G /, u' G I' in 2-adic camera images 1,1 = 
P 2 (Z 2 ), a 2-adic essential matrix E G Z 3 / 3 has been 
produced. This matrix is the limit of matrices E v G 
Z^? 3 . Through its factorisation 

E v = T V R V 

into a skew-symmetric matrix T v and a rotation i?„ it 
allows to determine a point U u G P 3 (Z 2 ^) of which the 
two cameras obtain the image points u Ul u' v mod 2". 
These are related through E v : 

u T v ■ E v ■ u' v = 0. 

The U v converge 2-adically to a point U G P 3 (Z 2 ), and 
the cameras 

k, k': P 3 (Z 2 ) -^F 2 (Z 2 ) 
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yield k(U) = u, k'(U) = u' . In the same way as u v 
is an approximation of pixel u at resolution v, the 3D- 
point U v is an approximation of voxel U at resolution v, 
because we can apply the map 12 in precisely the same 
manner also to three dimensions to obtain the inclusion 
P 3 (Z 2 )^ [0,2) 3 . 

5 Reconstruction from point correspondences 

We will in the following content ourselves with solving 
the reconstruction problem with two calibrated cameras 
by finding (possible candidates for) essential matrices 
in the 2-adic setting. 

Given n correspondences, we arrive at the system of 
linear equations 

ujEu'i, i = 1, . . . ,n. (3) 

For n < 9, we can use the Hensel lifting method to solve 
these linear equations, written as 

Ax = (4) 

with coefficient matrix A G Z^* 9 . Due to finite resolu- 
tion and our coding method, the matrix entries lie in 
some set {0, . . . , 2 V — 1}. Hence, the matrix does not 
change its shape modulo 2". 

Assume that we are given a basis B = (pi, . . . , b m ) 
in normal form of the solution space of (j4]) modulo 2. 
We will interpret B as an ordered set, but also as a 
matrix whose columns are the vectors b\, . . . ,b m . 

Theorem 2 (Linear Hensel's lemma) If the rank m 
of A mod 2 equals Rank(v4) ; then B lifts to a set of lin- 
early independent solutions of with natural numbers 
as entries. 

Proof The linear equation ([4]) is equivalent to a linear 
equation in staircase normal form. This normal form 
can be viewed as a system of Rank(A) equations in 
the 9 unknowns. By assumption, m — Rank(A). Hence, 
Theorem [TJ applies, and B lifts to a set of m solution 
vectors 

B = 1 61, . . .,b m \ 

to Q. 

Let us re-examine the proof of Theorem [T] in order 
to see that B is linearly independent. First, notice that 
B is read off the staircase normal form of A mod 2. 
This means that each b G B has in some row jb entry 
1, and below it all entries are zero. Further, all other 
entries of row jb of matrix B are zero. The sequence 
{jb)beB is strictly increasing with respect to the order 
of occurrence of b in B. Now, a lift of b to M fe ) modulo 



2 fc has the property that the entry in row jb is odd, and 
all entries below are even. Likewise, all other entries of 
row jb in a lift of B to B^ modulo 2 k are even. This 
description shows that the rows given by the sequence 
(jb) form a submatrix of B^ having rank m. Since 
this holds true for all k > 0, it follows that B is linearly 
independent. 

The lifting process stops after a finite number of 
steps, because of our initial assumptions. 

5.1 Reconstruction from 8 points 

Assume now that the rank of A mod 2 be n, and that a 
basis B — {61, . . . , bg- n } of the solution space be given 
in normal form. Then, by Theorem [2l finding a lift to 
Z2" of B yields a basis of the solution space of (j4|). 

As an application, we obtain a 2-adic version of the 
8-point algorithm of [TT], simply by setting n = 8. How- 
ever, in the same way as its original, this ignores the 
rank constraint det(-B) = for the essential matrix. 
Hence, we obtain the result: 

Theorem 3 Under the assumptions above for the ma- 
trix A G , the corresponding 8-point problem has 
a unique solution E, if additionally Rank(_B) — 2. If it 
can be assumed that A G Z 2 n , then E is computed after 
N — 1 iterations from E mod 2. 

Proof Since Rank(A mod 2) = 8, there is one solution 
basis vector b G Z 9 . By Theorem [21 it lifts to a non- 
trivial solution e of Ax = 0. If the the matrix E corre- 
sponding to e is of rank 2, it is the unique solution to 
the 8-point problem. The last assertion is an immediate 
consequence of the proof of Theorem [2] 

5.2 Reconstruction from 7 points 

The 7-point method by [6j[9] yields 7 linear constraints 
((3|) plus the cubic constraint 

det(E) = 0. (5) 

Let us write that system of equations as f(X) = 0. If 
the 7 points are sufficiently in general position, then the 
rank of Df (x) mod 2 is 8 for some solution x modulo 
2: 

f{x) = mod 2, (6) 
and we can lift to a 2-adic solution. The reason is: 
Lemma 1 For h := det(E) it holds true that 
D h (e)^0 mod 2, (7) 
if e G Z?, x3 is a sufficiently general instance of E. 
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Proof The polynomial h is of degree three in the vari- 
ables given by the entries of E. Let x be such an entry. 
Then h x := ^ is a polynomial of degree 2. Modulo 2 
only those terms of h x vanish which are of the form 2xy 
for some other variable j / i. Since h — det(E) is a 
square-free polynomial, this can never happen. 

In 2-adic analytic geometry, the conditions (|6|) and 
([7]) define an open subset in the space of all 3x3- 
matrices with entries in Z2. This can be seen easily by 
translating ([6]) to the inequality 

l/(*)l 2 < 1: 

and dealing similarly with J7])- 

We can go further and derive explicit conditions for 
the existence of ^-rational solutions. Namely, write 
down the 1-parameter solution of the linear equations 
as 

E = xE 1 + (l-x)E 2l 
and obtain a polynomial 

det(E) = h(x) = ax 3 + bx 2 + cx + d 

of degree 3. Now consider h(x) mod 2. In case d = 
mod 2, zero is a simple zero in F2 if and only if 

c = 1 mod 2. (8) 

In case d = 1 mod 2, one is a simple zero in F2 if and 
only if 

b = c = mod 2, (9) 
in case a = 1 mod 2, and 

6 = 0, c = 1 mod 2 (10) 
otherwise. 

Theorem 4 The 7 -point problem has a 1i 2 -rational so- 
lution for many choices of 7 point correspondences. Con- 
cretely, this is the case if (E$, 0), or U0\) hold true in 
their respective cases, together with the requirement that 
the rank of E be precisely 2. 

Proof Solve the 7-point problem modulo 2, and lift when- 
ever possible as discussed above. 

Remark 1 Due to non-linearity of the constraints, we 
cannot expect anymore to be able to lift solutions modulo 
2 to solutions which are defined over the natural num- 
bers. In other words, generally, the iteration will never 
stop unless an order of resolution is specified. Then, a 
2-adically approximate solution will be obtained. This 
is not different from the classical situation over the real 
numbers. 



5.3 Solving the 5-point non-linear equations 

Given five corresponding pairs of image points yields 
([3]) with n = 5. If the rank of the corresponding matrix 
is five, the general solution can be written as 

E = xE 1 +yE 2 + zE 3 + wE 4 . (11) 

Inserting this into the trace condition [3] 

2 ■ EE T E - Tr&ce(EE T ) ■ E = (12) 

yields 9 cubic equations in four variables. We follow the 
easily understandable method of elimination via hidden 
variables used by dOj, and obtain the linear equation 

C(z)-X = 0, 

after setting w = 1. Here, 

X = (x 3 ,y 3 , x 2 y, xy 2 , x 2 ,y 2 , xy, x, 1) 

is the vector of monomials, and the entries of C{z) are 
polynomials in z. Now, one seeks z such that 

det(C(») = 0. (13) 

The left hand side turns out to be a polynomial of de- 
gree 10. 

We will use the same lax formulation of our theorem 
as for the 7-point problem, but will be more specific in 
the proof: 

Theorem 5 There exists a l^i-rational solution to the 
5-point problem for many choices of 5 corresponding 
pairs of points. 

Proof Let 

10 

g(z) :=det(C(z)) =^0,^. 

i=0 

We consider two cases. 

First, assume oq = mod 2. In this case, g(0) = 
mod 2, and z = is a simple zero modulo 2 if and only 
if di = 1 mod 2. 

Secondly, if a\ = 1 mod 2, then /(l) = mod 2 
if and only if the number of odd coefficients in g(z) is 
even. Since 

f'(z) = agz s + a^z 6 + a^z A + a^z 2 + ai mod 2, 

1 is a simple zero modulo 2 if and only if in addition 
the number of odd coefficients in f'(z) is odd. 
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Remark 2 Hensel's lemma can be interpreted in this 
context as ap-adic stability result. Namely, assume that, 
due to correspondence error, a given choice of n points 
yields perturbed equations 

f{x) + e(x) = 0, 

where f{x) contain the "true" coefficients perturbed by 
some noise coefficients contained in e{x) with 2-adic 
maximum norm 

||e|| 2 =max{|ei| 2 \i = 1, . .. ,m} <2~ N , 

where m is the number of noise coefficients e^. This 
means that the noise coefficients satisfy the congruence 

ei=0 mod 2 N . 

Hence, the first N iterations of Hensel lifting will lead 
to identical approximations to the solution of the un- 
perturbed equations 

f(x) = 0. 

As the 2-adic encoding comes from interval subdivi- 
sions, this observation implies a greater stability in com- 
parison with the classical approach over the real num- 
bers. In particular, the existence of a liftable solution 
modulo 2 is not affected by perturbations e with ||e|| 2 < 
A. This is definitively in contrast to the situation over 
the real numbers, where small perturbations in coeffi- 
cients can drastically change the number of real solu- 
tions. That issue is addressed for the 5-point relative 
pose problem e.g. in JT]/. 

6 Conclusion 

An encoding scheme for image pixels through hierarchi- 
cal interval subdivision is proposed. This allows a 2-adic 
encoding of pixels driven by geometry. As an applica- 
tion to stereo vision, the 8-, 7- and 5-point equations are 
formulated with coefficients from the ring %2 of 2-adic 
integers. These polynomial equations are solved using 
some multivariate forms of Hensel's lemma. The essen- 
tial matrices are obtained in the form of sequences of 
matrices modulo powers of 2, corresponding to 2-adic 
approximations to the exact solutions in Z2, whenever 
these exist. Conditions on coefficients of equations mod- 
ulo 2 implying the existence of Z2-rational solutions are 
derived. One feature of the hierarchical encoding is that 
the number of iterations in solving the linear parts of 
the equations is logarithmic in the number of pixels. 
Also the precision in the solution of the non-linear equa- 
tions is directly related to the desired resolution in 3D. 



An immediate consequence of Hensel's lemma is that p- 
adically small perturbations of the equations do not af- 
fect the first approximations to their solution. Further, 
the existence of liftable solutions is not affected by rela- 
tively large perturbations. This indicates a greater com- 
putational benefit from the 2-adic approach compared 
to the classical approach using computational complex 
algebraic geometry before discarding non-real solutions. 

Acknowledgements Sven Wursthorn is thanked for the intro- 
duction into this fascinating topic, and Boris Jutzi for multiple 
fruitful discussions. 
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